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Abstract 

We have carried out a search for Main Belt Asteroids (MBAs) co-orbiting 
with the large MBA Vesta and the dwarf pla net Ceres. Through improving 



the search criteria used in IChristoul (j2000bl ) and numerical integrations of 



candidate coorbitals, we have identified approximately 51 (44) objects cur- 
rently in co-orbital libration with Ceres (Vesta). We show that these form 
part of a larger population of transient coorbitals; 129 (94) MBAs undergo 
episodes of co-orbital libration with Ceres (Vesta) within a 2 Myr interval 
centred on the present. The lifetime in the resonance is typically a few times 
~ 10^ yr but can exceed 2 x 10^ yr. The variational properties of the or- 
bits of several co-orbitals were examined. It was found that their present 
states with respect to the secondary are well determined but knowledge of 
it is lost typically after ~ 2 x 10^ years. Objects initially deeper into the 
co orbital region mainta in their coorbital state for longer. Using the model 



of iNamouni et al.l (119991 ) we show that their dynamics are similar to those of 
temporary coorbital NEAs of the Earth and Venus. As in that case, the life- 
time of resonant libration is dictated by planetary secular perturbations, the 
inherent chaoticity of the orbits and close encounters with massive objects 
other than the secondary. In particular we present evidence that, while in 
the coorbital state, close encounters with the secondary are generally avoided 
and that Ceres affects the stability of tadpole librators of Vesta. Finally we 
demonstrate the existence of Quasi-satellite orbiters of both Ceres and Vesta 
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and conclude that decametre-sized objects detected in the vicinity of Vesta 
by the DAWN mission may, in fact, belong to this dynamical class rather 
than be bona-fide (i.e. keplerian) satellites of Vesta. 
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The coorbital resonance, where the gravitational interaction between two 
bodies with nearly the same orbital energy leads to stable and predictable mo- 
tion, is ubiquitous in the solar system. Objects attended by known co-orbital 
companions include Jupiter, Mars a s well as the satu rnian satellites Tethys, 
Dione, Janus and Epimetheus (see IChristoul . l2000a, for a review). More 



recently, the planet Neptune was added to this list 



Sheppard and Trujillo 



20061 ) w hile an additional co orbital of Dione was discovered by the Cassini 



mission (IMurray et al.l . 120051 ). In all these cases, the motion h as been shown 



to be stable against all but the most slow-acting perturbations (ILissauer et al. 



1985: Levison et al.. 1997: Brasser et al.. 2004 



Scholl et al.l . l2005[ ). 



The discovery of a coorbital att e ndant of th e Earth on a highly inclined 
and eccentric orbit ( Wiegert et al. . Il997l Il998l ) motivated new theoretical 
work in the field. NamounT l 19991 ) . using Hill's approximation to the Re- 
stricted Three Body Problem (R3BP) showed analytically that the introduc- 
tion of large eccentricity and inclination modifies considerably the topology of 
coorbital dynamics near the secondary ma ss. It results in the appearan ce of 
bounded eccentric orbits (Quasi-Satellites; iMikkola and Innanenl . 119971 ) and, 
in three dimensions, "compound" orbits and stable transitions between the 
different modes of libration. Further, he demonstrated numerically that these 
results hold when the full R3BP is considered. In this case, the appearance 
of new types of compound orbits, such as asymmetri c modes or com p ounds 
of tadpoles and retrograde satellites, was shown in iNamouni et al.l (119991 ) 



to be due to the secular evolution of the coorbital potential. Such types of 
coorbital libra t ion w ere identified in the motion of the object highlighted by 



Wiegert et al.l fll997l ) as well as other near-Earth asteroids but the secular 



forcing o f the potential in that case is provided by planetary secular pertur- 



bations (INamouni et al.l . Il999l : IChristoul . l2000al ). The expected characteris 



tics of the population of co-o rbita ls of Earth and Venus were in vestigated by 



Morals and Morbidellil (120021 ) and lMorais and Morbidellil (120061 ) respe ctively. 



Christoul (l2000bl ). motivated by Pluto's ability, as demonstrated by lYu and Tremaine 
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(Il999l ) and iNesvorny et al.l (120001 ) . to trap other Edgeworth-Kuiper Belt ob- 
jects in co-orbital motion with itself, demonstrated in turn that main belt 
asteroids can co-orbit with the dwarf planet 1 Ceres and the large Main Belt 
Asteroid (MBA) 4 Vesta. Four such asteroids were identified, two co-orbiting 
with Ceres and two with Vesta. 

Here we report the results of a search for additional coorbitals of these 
two massive asteroids. This was motivated partly by the large growth in the 
number of sufficiently well-known MBA orbits during the intervening decade, 
but also a refinement of the search criterion used in the work by Christou. As 
a result, we find over 200 new transient co-orbital MB As of Ceres and Vesta. 
In this work we examine their ensemble properties, use existing dynamical 
models to understand how they arise and identify similarities with co-orbital 
populations elsewhere in the solar system, in particular the transient coor- 
bital NEAs of Earth and Venus. 

The paper is organised as follows: In the next Section we expose our 
Search method ology, in p a rticula r those aspects which differ from the search 
carried out by IChristoul ( l2000bl ). In Section 3 we describe the statistics 
of coorbital lifetime and orbital element distribution found in our integra- 
tions. In addition, we examine the robustness of the dynamical structures 
w e observe. In Section 4 we investigate the degree to which the model of 



Namouni et al.l (119991 ) can reproduce the observed dynamics. Section 5 deals 



with the effects of additional massive asteroids in the three-body dynamics 
while Section 6 focuses on the stability of so-called Quasi-Satellite orbits. 
Finally, Section 7 summarises our conclusions and identifies further avenues 
of investigation. 



58 



60 



2. Search Methodology 



Christou! (l2000bl ) searched for candidate Ceres coorbitals by employing 



59 the osculating semimajor axis of the asteroid relative to that of Ceres 



(equal to (a — accrcs) /accrcs) to highlight objects that merited further in- 

61 vestigation. This method, although it led to the successful identification of 

62 two coorbitals, 1372 Haremari and 8877 Rentaro, ignores the existence of 

63 high frequency variations in a due to perturbations by the planets, especially 

64 Jupiter. This is illustrated in the upper panel of Fig. [T] where the evolution 

65 of flr for asteroid 1372 Haremari, one of the objects identified by Christou, is 

66 depicted. A periodic term of amplitude 0.001 apparently causes to move 

67 in and out of Ceres' coorbital zone (bounded by the dashed horizontal lines) 
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68 every 20 years. On the other hand, an object with < e would have a con- 

69 junction with Ceres every > 6 x 10^ years. Here e = denotes Hill's 

70 parameter in the Circular Restricted Three-Body Problem for a secondary 

71 body with a mass /i scaled to that of the central body (in this case, the Sun). 

72 This suggests that coorbital motion should be insensitive to these variations 

73 and that a low-pass filtered semimajor axis would work better as an indicator 

74 of libration in the 1:1 resonance with Ceres. We have chosen this to be the 

75 synthetic proper ser nimajor axis (hereafte r refer red to as "proper semimajor 
axis") as defined by lKnezevic and Milani fl200(t ). Synthetic proper elements 



76 



100 



77 are numerically-derived constants of an asteroid's motion under planetary 

78 perturbations. The proper semimajor axis, in particular, is invariant with 

79 respect to secular perturbations up to the second order in the masses under 

80 Poisson's theorem and, hence, an appropriate metric to use in this work. In 

81 the bottom panel of Fig. [1], we plotted the proper relative semimajor axis 

82 ttr^ p = [ttp — Ceres) /op, Ceres; bold horizontal line) superimposed on the 

83 ttr history of the asteroid but filtered with a 256-point running-box average. 

84 The slopes at the beginning and the end of this time series are the result 

85 of incomplete averaging. The dots in the bottom panel represent the un- 

86 averaged osculating semimajor axis sampled with a constant timestep of 1 

87 yr. It is observed that the proper semimajor axis agrees with the average of 

88 its osculating counterpart, whereas the osculating value sampled at t = ( 

89 Julian Date 2451545.0; dashed horizontal line) does not. 

90 For our search we used the database of synthetic proper elements for 

91 185546 numbered asteroids computed by Novakovic, Knezevic and Milani 

92 that was available as of 09/2008 at the AstDys online information service 

93 (http:/ /hamilton. dm. unipi.it/~astdys/propsynth/numb.syn). A total of 648 

94 and 514 main-belt asteroids respectively with proper semimajor axes within 

95 ±ae of those of Ceres and Vesta were identified. Table [T] shows the pa- 

96 rameters relevant to this search for the massive asteroids. The second row 

97 provides the mass ratio fi of the asteroid relative to the mass of the Sun, the 

98 third row its Hill parameter e, the fourth row the asteroids' synthetic proper 

99 semimaxor axis as given in the database and the fifth row the product 
, X e which is the half-width of the coorbital region in AU. 



101 The state vectors of these asteroids at Julian Date 2451545.0 were re- 



102 trieved from HORIZONS flGiorgini et al.l . Il996l ) and numerically integrated 



103 one million years in the past and in the future using a model of the solar 

104 system consisting of the 8 major planets and the asteroids 1 Ceres, 4 Vesta, 

105 2 Pallas and 10 Hygiea. Mass values adopted for these four asteroids were 
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taken from iKonopliv et al.l (j2006|). Their initial state vectors as well as plan- 
etary initial conditions and constants were also retrieved from HORIZONS. 
The integrations were carrie d out using the " hybrid" scheme which is part 
of the MERCURY package fIChambersl . Il999h . This scheme is based on a 
second-order mixed variable symplectic (MVS) algorithm; it switches to a 
Bulirsch-Stoer scheme within a certain distance from a massive object. For 
all the integrations reported here, this distance was 5 Hill radii. A time step 
of 4 days or ~ l/20th of the orbital period of Mercury was chosen in orde r 
to mitigate the effects of stepsize resonances (IWisdom and Holm ani. Il992h. 
Trials of this scheme vs an RA15 RADAU integrator (lEverhart . 1985) in- 
cluded in the same package and with a tolerance parameter of 10"^'^ showed 
the results to be indistinguishable from each other while the hybrid scheme 
was significantly faster. 



119 3. Results 

120 3.1. Population Statistics 

121 In our runs, we observed a total of 129 and 94 asteroids enter in one of 

122 the known libration modes of the 1:1 resonance with Ceres or Vesta respec- 

123 tively at some point during the integrations. A full list of these asteroids 

124 is available from the corresponding author upon request. Table [2] shows a 

125 statistical breakdown of the observed population according to different types 

126 of behaviour. The first row identifies the secondary (Ceres or Vesta). The 

127 top part of the Table shows the number of asteroids that were captured into 

128 co-orbital libration at some point during the 2 x 10^ yr simulation (second 

129 column), the number of asteroids that were found to be co-orbiting at present 

130 ( "current" coorbitals; third column) and the number of current L4 and L5 

131 tadpole librators (fourth and fifth rows respectively). The bottom part of 

132 the table shows the number of current horseshoe librators (second column), 

133 the number of objects currently in transition between two distinct libration 

134 modes (third column) and the number of current Quasi-Satellite (QS) libra- 

135 tors. The fifth column shows the number of asteroids that were librating 

136 for the full duration of the numerical integration, although not necessarily 

137 in the same mode ("persistent" coorbitals). Numbers in brackets refer to 

138 MBAs that remained in the same libration mode throughout the integration 

139 ("single mode persistent" coorbitals). A total of 95 MBAs are currently in 

140 the 1:1 resonance with either Ceres or Vesta with the ratio of objects cor- 

141 responding to the two asteroids (51/44) being roughly the same as for the 
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142 original sample (129/94). The number of MBAs currently in tadpole orbits 

143 is 19 and 12 for the two asteroids respectively while horseshoes are evenly 

144 matched at 20 and 21. 11 objects for each asteroid are in the process of tran- 

145 sitioning between different modes of libration. One object, MBA 76146, is 

146 currently in a Quasi-Satellite (QS) orbit around Ceres with a guiding centre 

147 amplitude of ~ 10° in (rightmost panels in Fig [2]) while two other MBAs, 
138245 and 166529, are currently transitioning into such a state. 



148 



150 



3.2. Some Examples 

In most cases, the librations were transient i.e. changing to another mode 

151 or to circulation of the critical argument = A — Aceres/Vesta- In that sense, 

152 there is a similar i ty to what is observed for co-orbitals of the Earth and Venu s 
fiNamouni et all . Il999l : IChristoul . boOOat iMorais and Morbidelli l2002l . l2006h . 



153 



154 Compound modes (unions of Quasi-Satellite and Horseshoe or Tadpole li- 

155 brations) also appear, although they are rare. Examples of different types of 

156 behaviour are shown in Figs [2] and E] for Ceres and Vesta respectively. 

157 In the former case, asteroid 71210 (left column) is currently (t = at Ju- 

158 lian Date 2451545.0) in a horseshoe configuration with Ceres, transitioning 

159 to tadpole libration at times. The position of Ceres at A^ = is avoided. 

160 Asteroid 81522 (centre column) librates around the L4 triangular equilibrium 

161 point of Ceres for the duration of the integration. In the right column we 

162 observe stable transitions of the orbit of asteroid 76146 between different 

163 libration modes. The maximum excursion of the guiding centre of the as- 

164 teroid from that of Ceres during the integration is ar = 10~^ or 0.2e, well 

165 within the coorbital region. This asteroid becomes a Quasi-Satellite of Ceres 

166 at t = — 2 X 10^ yr and transitions into a horseshoe mode at t = +10^ yr. 

167 In the case of Vesta, asteroid 22668 (left column) transitions from a pass- 

168 ing to a horseshoe mode and back again. It is currently a horseshoe of Vesta. 

169 Near the end of the integration it enters into Quasi-Satellite libration with 

170 Vesta, where it remains. At t = +5 x 10^ yr it executes half a libration in a 

171 compound Horseshoe/ Quasi Satellite mode. Asteroid 98231 (centre column) 

172 is an L5 tadpole of Vesta at t = — 1 x 10^ yr. The libration amplitude increases 

173 gradually until transition into horseshoe libration occurs at t = —3 x 10^ yr. 

174 During this period the libration amplitude begins to decrease until the re- 

175 verse transition back into L4 libration takes place at t = +6 x 10^ yr. This 

176 behaviour should be compared with the case of asteroid 71210 in Fig. [2l The 

177 evolution of asteroid 156810 (right column) is similar to that of 98231 except 
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178 that the rate of increase of the hbration amphtude reverses sign before tran- 

179 sition into horseshoe hbration takes place. As a result, the asteroid librates 

180 around the L4 equilibrium point of Vesta for the duration of the integration. 

181 3.3. Variational Properties 

182 Although these results can be regarded statistically, their value is in- 

183 creased if we establish their robustness against the ephemeris uncertainties 

184 of the asteroids. To investigate this, we have picked 6 asteroids in the sam- 

185 pie, populated their 1-sigma uncerainty elhpsoids with 100 clones per object 

186 using states and state covariance matrices retrieved from the AstDys online 

187 orbital information service, and propagated them forward in time for 10^ yr 

188 under the same model of the solar system as before but using an integrator 

189 developed by one of us (PW) with a 4-day time step. This code utilises the 

190 same algorithm as the HYBRID code within MERCURY. In the original inte- 
rn grations, 65313 and 129109 persist as L5 Trojans of Ceres, 81522 and 185105 

192 as L4 Trojans of Ceres and 156810 persists as an L5 Trojan of Vesta. Asteroid 

193 76146 is a Quasi-Satellite of Ceres until -|-10^ yr. Their Lyapounov times, 

194 as given by the proper element database of Novakovic et al, are: 3 x 10^, 

195 5 X 10^, 2 X 10^, 6 X 10^, 10^ and 5 x 10^ yr respectively. 

196 The results of this exercise are shown in Fig. |H We find that all clones of 

197 65313 and 81522 remain as Trojans of Ceres for the full integration. 76146 

198 persists as a Quasi-Satellite of Ceres until ~ -|-6 x 10^ yr. Three clones of 

199 129109 escape from hbration around L5 after ~ 5 x 10^ yr, but most remain 

200 in hbration until the end. Those of 185105 begin to diverge at ~ -|-3 x 10^ 

201 yr until knowledge of the object's state is lost near the middle of the in- 

202 tegration timespan. The clones of the Vesta Trojan 156810 suffer a simi- 

203 lar fate, but divergence is more general in nature; most of the clones enter 

204 into other hbration modes by -|-6 x 10^ yr. We conclude that the present 

205 state of these objects, as determined by the original integrations, is robust. 

206 The differences we observe in their evolution could be due to several causes. 

207 Firstly, although we use the same force model and the same type of inte- 

208 grator, the initial states of the asteroids in the new integrations were taken 

209 from AstDys, not HORIZONS, and correspond to a later epoch of oscula- 

210 tion. In addition, the volumes of space sampled by the object's ephemeris 

211 uncertainties, being proportional to the eigenvalue product of the respective 

212 covariance matrices, are smallest for 65313 and 81522. We therefore expect 

213 that, assuming similar Lyapounov times, clones of low- numbered asteroids - 

214 generally those with longer observational arcs - would be slower to disperse 
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215 than those of high-numbered ones. Interestingly, the clones disperse faster 

216 than one would expect from their Lyapounov times, however some differences 

217 would be expected since the model used to compute these did not include the 

218 gravitational attraction of large asteroids. In Section 5, the role of impulsive 

219 perturbations during close encounters is investigated in detail. 

220 3.4- Proper element distribution 

221 In order to understand the ensemble properties of the population and 

222 how these might differ from those of other MBAs, we have compared the dis- 

223 tribution of their proper elements to that of the broader population. Fig [5] 

224 shows the distribution of the relative proper semimajor axis (blue curve) 

225 ttr^ p = {ap- Ceres/Vesta) / Op, Ceres/Vesta of thcSC MBAs, SCaled tO €, SUper- 

226 posed on that of all MBAs in the interval [— 2e, +2e] (red curve). The width 

227 of each bin for both panels is 0.06. The sharp peak at = —1.2 in the 
22B plot for Vesta reflects an increased concentration of asteroid proper elements 

229 probably associated with the 36:11 mean motion resonance with Jupiter. Al- 

230 though interesting in its own right, we do not, at present, have reason to 

231 believe that its existence near the co-orbital region of Vesta is anything more 

232 than coincidental. Hence, we refrain from dicussing it further in this paper. 

233 In both cases, the distributions appear to be centred at a^, p = 0. Gaus- 

234 sian fits to the centre fi and the standard deviation a of the distribution give 

235 a Full Width at Half Maximum (FWHM; 2^2 \og2a) of 0.334 ± 0.033 and a 

236 centre at —0.053 ib 0.017 for the Ceres distribution. The slight offset to the 

237 left is probably due to slightly higher counts for the bins left of a^, p = 0. 

238 The Vesta distribution is slightly narrower (FWHM of 0.292 ± 0.028) but 

239 more symmetric around the origin (centre at —0.006 ± 0.014). No cases of 

240 coorbital libration were observed for asteroids with ja^, p| > 0.42 while only 

241 two cases (both with Vesta) had la^, p| > 0.30. 

242 Seeking additional insight into the dependence of the coorbital state on 

243 the semimajor axis, we examined the distributions of different types of co- 

244 orbitals - as observed in our simulations - normalised to the total number 

245 of MBAs in each bin. In Fig. [6] we show the distributions of all coorbitals 

246 (red curve), current coorbitals (blue curve) and persistent coorbitals (gray 

247 curve). Fitted values of the Gaussian parameters (/i, FWHM) for the three 

248 populations are given in Table [31 The distributions for Vesta coorbitals 

249 are consistently narrower than those of Ceres implying that this is a real 

250 difference between the two. 



8 



251 Interestingly, both cases exhibit a hierarchy in the three populations: the 

252 persistent population is embedded into the current one, which in turn is 

253 embedded into the distribution of all objects. One important consequence 

254 of this observation is that one can robustly define the boundaries of each 

255 population. This is, of course, partly due to the criteria used to define each 

256 population but the fact that there are clear differences between the three 

257 populations (ie no two populations coincide) is not a trivial one. Hence, 
persistent coorbitals are confined in the domain |ar, p\ < 0.12, current coor- 



258 



259 bitals within ja^, p\ < 0.24 and all coorbitals in the domain \ar^ p\ < 0.42. 



260 



The shape of the distributions in Figs. [5] and [6] can be partly attributed 

261 to the coorbital dynamics (see Section 4). However, they must also be af- 

262 fected by chaotic 'nois e' in the determination of t he proper elements which 



279 
280 



284 
285 
286 
287 
288 



264 tions, 



263 we us ed in our search ( Milani and Knezevid . Il994l ) . In their 2 Myr Integra- 



Knezevic and Milanil ( 120001 ) regarded the derived proper elements of 

265 MBAs with dap < 3 X 10"^, cTg^ < 3 x 10"^ and ai^ < 10'^ as 'good'. All but 

266 four of the asteroids considered here belong to this catelogy. The bound for 

267 the proper semimajor axis corresponds to 20% of the width of the coorbital 

268 region of Ceres and 35% of that of Vesta (Table [T]). It is also comparable 

269 to the fitted widths of the distributions of the coorbitals found here. Hence, 

270 the actual distributions are likely significantly altered by a convolution with 

271 an error function. On the other hand, this convolution does not completely 

272 smear out the true distribution of a^, p since, in that case, the observed sort- 

273 ing of the populations according to residence time in the resonance would 

274 not occur (Fig. [6]). 

275 The distribution of the proper eccentricities and inclinations of individual 

276 coorbitals in relation to those of other MBAs are shown in Fig. [71 Plus sym- 

277 bols denote MBAs that have tested negative for co-orbital motion within the 

278 period [ — 10® yr, 10® yr] . Asterisks and squares refer to the respective pop- 
ulations of current and persistent coorbitals while the filled circle marks the 
location of the secondary (either Ceres or Vesta). In the interests of clarity, 

281 we have not plotted the distribution of all coorbitals. Instead, we show as 

282 triangles those persistent coorbitals that remained in either L4 or L5 tadpole 

283 libration for the full simulation. Their location deep into the coorbital region 
are in agreement with the theoretical upper limit - ^y (8/3)yUccrcs/Vcsta - for 
near-planar, near-circular tadpole orbits which evaluates to 0.065e for Ceres 
and 0.053e for Vesta. 

The significant size of the sample of coorbitals under study prompted a 
search for trends in the distribution of A^. Fig. IHl shows a histogram of this 
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289 quantity at t = for all current coorbitals and with a bin size of 30°. On this 

290 we superimpose, as a dashed line, a histogram of all objects which are not 

291 currently coorbiting with either Ceres or Vesta. The vertical line segments 

292 indicate square-root Poissonian uncertainties. There appear to be no features 

293 that stand out above the uncertainties. Hence the coorbital resonance does 

294 not measurably affect the phasing of the populations of current coorbitals 

295 with respect to their secondary in this case. 



296 4. Analysis of the Dynamics 



297 



The dynamical context presents some similarities with coorbitals of the 

298 Earth and Venus such as non-negligible eccentricities and inclinations. Here 

299 we attempt to model the evolution using the framework of the restricted 

300 three body problem where a p article's state evolves under the gravity of the 



301 

302 



304 



307 
308 



310 



312 



Sun and the secondary mass (INamouni et al.l . Il999l ). This is done through 
the expression 

al = C-^S{K,eJ,u) (1) 



303 where 



s=^ r (-^-'-^]dx. (2) 

271 J^^ V|r-rs| aas ' ^ ^ 



Here, a, e, /, u and A denote the semimajor axis, eccentricity, inclination, 

305 argument of pericentre and mean longitude of the particle's orbit. The sub- 

306 script "S" is used to denote the same quantities for the secondary. The 
heliocentric position vectors of the particle and the secondary are denoted as 
r and rg respectively. The relative elements and are defined as 



ar = {a - as) /as, Xr = >^ - (3) 

309 and fi is the mass of the secondary scaled to the total system mass. Eq. [1] 
can be seen as a conservation law where C is the constant "energy" of the 



311 particle, af its k inetic energy and the term containing S its potential energy. 



Namouni et al.l showed that, as the left-hand side of this expression cannot 

313 be negative, it restricts, in general, the evolution of {ar,Xr)- A collision 

314 (r = rs) can only occur for specific combinations of values for e, / and u. 

315 Hence, actual collisions are rare and the above formulation is generally valid. 

316 For computatio nal purposes, Eq.[2]may be e valuated using standard two- 



317 body formulae (eg [Murray and Dermottl . Il999l ) as A 3> e, /, cu and A^. One 
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318 other consideration that is specific to this paper is that the high frequency 

319 harmonics of a are external to the R3BP (hence external to the model) and 

320 must somehow be removed before the above expressions may be used. How- 

321 ever, at t = 0, the relative proper semimajor axis a^, p can be considered to 

322 be this low-pass filtered value of ar in the Sun - Ceres/ Vesta - MBA problem 

323 (see also Section 2). Hence, the constant C can be evaluated and the model 

324 can be readily applied. In Fig. IHlwe show S profiles for each of the MBAs in 

325 Figs [2] and [3] compared to the location of the point (A^, £ = 3C/(8/i)) at t = 

326 (dotted circle). As motion is restricted to the domain above S, the model 

327 predicts that the Ceres coorbitals 71210, 81522, 76146 are currently in L5 

328 tadpole, L4 tadpole and Quasi-Satellite libration respectively. Similarly the 

329 Vesta coorbitals 22668, 98231 and 156810 are predicted to be in horseshoe, 

330 horseshoe and L5 libration respectively. Referring to the Figure, the model 

331 apparently succeeds in 5 out of the 6 cases, but fails in the case of the Ceres 

332 coorbital 71210 where the observed mode of libration is a horseshoe. 

333 This is probably due to the fact tha t Eg [1] is evalua t ed wh en ar = i.e. 



334 at the turning points of the libration. In lNamouni et al.l (119991 ). ^ e, I, u 

335 and the orbit can be considered "frozen" during a libration cycle. In our case, 

336 however, we observe that cj > A,, because of the small mass of the secondary. 

337 Incidentally, this parity in timescales may also account for the general lack 

338 of compound libration modes for these coorbitals, as u controls the relative 

339 height of the maxima of S on either side of A^ = 0. 

340 To quantify the effect that this has to our model, we evaluated S against 

341 S for the example MBAs shown in Figs [2] and [3] but at different values of Ar 

342 and u. We found that determination of S is generally insensitive to u, except 

343 near the local maxima bracketing the origin on the A^ axis. Physically, these 

344 correspond to the closest possible cartesian distances between the particle 

345 and the secondary so it is not a surprise that they are sensitive to the orbital 

346 elements. Particularly for the case where the model failed, £^ — S" ~ 0.6 when 

347 the asteroid reaches the far end of the model tadpole (A^ ~ —130°) and the 

348 potential maximum at Ar = 180° i.e. the object is classified as a horseshoe in 

349 agreement with the numerical integrations. Hence, this method for determin- 

350 ing the resonant mode is formally valid where the object is currently near the 

351 turning point of the libration i.e. those of MBAs 81522, 76146 and 156810. 

352 In the other cases, the more involving process of monitoring the quantity 

353 S — S in the integrations for a time period comparable to a libration cycle 

354 would be necessary to establish the libration mode. 

355 Finally, we wish to understand the stability of the QS librator of Ceres, 
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356 76146, in the context of our findings. The sensitivity of 5" to w for the 

357 local maxima near = is important as these features of the potential 

358 are the "gatekeepers" for evolution in and out of the QS. This object was 

359 captured as a QS from a passing orbit and completes 7 cycles in this mode 

360 before becoming a horseshoe. It is currently at — 8°. Irrespective of the 

361 value of escape from the QS is certain when S't=o < 'S'a^=o i-e. when the 

362 extremum at the origin becomes a local maximum. According to the model 

363 this occurs for —34° < a; < +12°. Keeping in mind the vr periodicity of S in 

364 w, the asteroid will turn back 360/2/46 ~ 4 times before it escapes, in good 

365 agreement with what is observed in the numerical simulation. 



366 5. The role of other massive asteroids 



367 



381 
382 
383 



388 
389 



Christoul (j2000af ) showed that Venus and Mars play a key role in the 

368 evolution of Earth co-orbitals. These can force transitions between different 

369 libration modes or escape from the resonance altogether. Here the only 

370 candidates available to play a similar role are other massive asteroids. In 

371 this paper we have focused on the effects of Ceres and Vesta - as well as 

372 Pallas - on Vesta or Ceres co-orbitals respectively. In a first experiment to 

373 determine their role (if any), we have integrated the same six asteroids as 

374 in Section 3 but, in the first instance leaving Pallas out of the model ("No 

375 Pallas" or NP) and, in the second, only under the gravity of the secondary (ie 

376 Ceres or Vesta as appropriate; "Secondary Only" or SO). We find that the 

377 evolution of the Ceres Trojans 65313 and 81522 and of the Quasi-Satellite 

378 76146 are the same in both of these runs as well as the original runs where all 

379 three massive asteroids were present ("All Masses" or AM). In contrast, we 

380 find significant differences in the evolution of the remaining three asteroids, 
129109, 156810 and 185105 (Fig.flO]). In the SO runs (top row), all clones of 
these Trojan librators persist as such. This is not the case for the NP and AM 
runs (middle and bottom rows respectively). Most clones ultimately leave 

384 this libration mode, with the exception of 129109 in the AM runs. Although 

385 it is difficult, on the basis of a sample of six, to attempt to decouple the 

386 effects of the individual massive asteroids, these results indicate that their 

387 presence does have an effect on the lifetime of the co-orbital configurations 
we have found. 

To quantify these in a more statistically robust sense, we used MERCURY 

390 as before to re-integrate the original 129 and 94 co-orbital MB As of Ceres 

391 and Vesta respectively under the SO, NP and AM models for 10^ yr in the 
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392 future (ie half the timespan of the original integrations). We then check which 

393 MBAs survive, either as "persistent" or "single mode" libration in the new 

394 runs. The results are shown in Table |H As far as persistency is concerned, 

395 the fraction of those co-orbitals in the respective samples is between 20% and 

396 25% for the three models examined (SO, NP and AM) but also for both Ceres 

397 and Vesta co-orbitals. Upon closer inspection, we find that approximately 

398 the same number of persistent co-orbitals (16 for Ceres and 14 for Vesta) 

399 are common between the three models. Thus, Vesta harbours a slightly 

400 larger fraction of these co-orbitals than does Ceres. The remainder of the 

401 persistent population is composed of objects that gain or lose the persistency 

402 property from one model to the next. There is no apparent trend towards one 

403 direction or the other; objects that lose this property are replenished by the 

404 objects which gain it. Thus, it seems that persistency of co-orbital motion is 

405 independent of the adding or removing massive asteroids from the model. If 

406 this premise is correct, then a mechanism that can explain the observations 

407 is the intrincic chaoticity of the orbits. We saw in Section 3 that the amount 

408 of chaotic 'noise' in the proper semimajor axis is comparable to the width of 

409 the coorbital zones of Ceres and Vesta. It is thus reasonable to expect that 

410 some objects are removed from the resonance while others are injected into 

411 it as a consequence of that element's random walk. 

412 The behaviour of single- mode coorbital libration is altogether different. 

413 For Ceres, we find that the addition of other massive asteroids increases the 

414 number of persistent single-mode librators, from 6 (SO) to 14 (NP, AM). For 

415 Vesta, we observe the opposite trend: adding Pallas and Ceres to the model 

416 decreases the number of such librators from 11 (SO) to 3 (NP) and 5 (AM). 

417 In addition, only one object, 87955, persisted as a Ceres L4 Trojan in all 

418 three models. The statistics are marginally significant for these low counts. 

419 Nevertheless, they compelled us to explore possible causes. 

420 It is tempting at this point to try and correlate persistency with the 

421 objects' Lyapounov Characteristic Exponent (LCE), an established quantifier 

422 of chaoticity, available from the proper element database used in this paper. 

423 However, we point out that these LCEs were computed under a dynamical 

424 model that does not involve any massive asteroids. As it is not clear how 

425 the presence of the co-orbital resonance will affect the determination of the 

426 LCE, we refrain from attempting such a correlation in this work. 

427 In Fig. [11] and [12] we present the distribution of the number of Ceres 

428 or Vesta coorbitals respectively that undergo a given number of encounters 

429 within 5 Hill radii {Rh = ae) of the massive asteroids in the three models. 
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A number of interesting features are evident. Firstly, persistent coorbitals 
(green and blue columns), as well as a certain fraction of non-persistent coor- 
bitals (red columns; 25-30% for Ceres coorbitals, ~50% of Vesta coorbitals), 
d o not, in general, approach the secondary. This feature was also identified 



m 



Christoul (l2000al ) for Earth co-orbitals (cf Fig. 8 of that paper). It seems to 



be a generic feature of coorbital dynamics at high eccentricity and inclination 
orbits although the cause-and-effect relationship is not yet clear. 

In addition, the distribution of encounters for non-persistent coorbitals 
(red) shows a tail; in other words, many of these objects undergo many en- 
counters with the secondary. This is probably due to the slowness of the evo- 
lution of the relative longitude of the guiding centre compared to the epicyclic 
motion. Indeed, we find that these encounters are not randomly distributed 
in time but occur in groups typically spanning a few centuries. In contrast, 
those related to massive asteroids other than the secondary have an upper 
cutoff (5 encounters for Pallas and Vesta, 10 for Ceres). The shapes of the 
distributions for the three classes of coorbitals are similar. The apparently 
high frequencies observed in the case of Ceres encounters for single-mode 
persistent coorbitals of Vesta are probably due to the small size of the pop- 
ulation in those classes. The better populated distributions in the AM and 
SO models mimic the distributions of the other two classes we investigated. 

A different way to look at the data is to create histograms of minimum 
distances, since distance is one of the factors (velocity being the other) that 
determine the magnitude of the change in an MBA's orbit. Fig. [T3] and 
[Ml show the distribution of these distances in units of -Rh- Bin i contains 
all recorded encounter distances between [i — 1) and iRn and have been 
normalised with respect to the area of the corresponding annulus of width i?H 
on the impact plane. For Ceres coorbitals we do not discern any statistically 
significant variations, in other words the counts within the different bins are 
the same given the uncertainties. For Vesta coorbitals, the situation is similar 
with one exception: we note that single mode persistent coorbitals (blue) do 
not approach Ceres closer than one Hill radius. Due to the low counts, we 
cannot exclude the possibility that we are looking at statistical variation 
in the data. However, if it is the signature of a real trend, it would mean 
that encounters with Ceres may cause Vesta coorbitals to exit a particular 
libration mode. 

Regardless of their significance at population level, there is clear evidence 
in the data that Ceres encounters can affect the orbital evolution of coorbitals 
of Vesta. Fig [T^ shows two instances where Vesta coorbital MBAs 45364 and 
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468 164791 leave the resonance following close encounters with Ceres deep within 

469 the Hill sphere of that dwarf planet. We have searched for similar occurences 

470 in the evolution of Ceres coorbitals but without success. Neither Pallas nor 

471 Vesta appear capable of playing a similar role. This is probably due to their 

472 smaller masses (hence physically smaller Hill spheres) but also, specifically 

473 for the case of Pallas, the higher encounter velocities. 



6. Quasi-Satellites 



476 
476 
477 
478 
479 
480 
481 
482 
483 
484 
485 
486 
487 



This Section is devoted to the existence, as well as stability, of so-called 
"Quasi-Satellite" (QS) or "bound" orbits. These appeared first in the lit- 
erature as Retro grade Satelhte (RS) orbits (as it turned out, a special case 
of the QS state; Jackson . 19131) and later studied in the con text of dynami- 
cal systems analysis (Henon, 19691 : Henon and Guyot . 19701 ). More recently, 
the survival of this libra t ion rn ode i n the real sola r syste m was examined 
by iMikkola and InnanenI fll997f ) and IWiegert et aD fl2000h . Currently, two 
known Quasi-Satellit es exist for the Earth (Wajer, 20101 ) and one for Venus 
flMikkola et al.l . l2004h . 

Several instances of QS libration relative to Ceres and Vesta were found 
in our simulations. Before these are discussed in detail, we introduce some 
elements of a co nvenient t h eoret ical framework to study QS motion dynamics, 
namely that by iNamounil (119991 . see also Henon and Petit, 1986). It employs 
the set of relative variables 



Xr = Cr COS {ut — ZUr) + Ctr 
Ur = —2er sin {nt — ZUr) + Xr 

Zr = Ir sin {nt — Vt^) 



(4) 
(5) 
(6) 



489 where n is the mean motion and e^, /r, '^r and VL^ the relative eccentricity, 

490 relative inclination, relative longit ude of per i centr e and relative longitude of 

491 the ascending node as defined in iNamounil (119991 ) respectively. In this for- 

492 mulation, the motion is composed of the slow evolution of the guiding centre 

493 (or, \r) to which a fast, three-dimensional epicyclic motion of frequency n 

494 and amplitude proportional to and is superposed. Examples of guiding 

495 centre libration while in QS mode have been shown in Figures[2] (bottom right 

496 panel) and [3] (bottom middle panel). To illustrate the relationship between 

497 the two components of the motion, we show in Fig [16] two examples of QS 

498 motion recovered from our simulations. The left panel shows the cartesian 
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502 
603 



607 
608 
609 



499 motion of MBA 50251 in a cartesian heliocentric frame that rotates with the 

600 mean motion of Ceres ('C') and for a period of ~ 5 x 10^ yr. The guiding 

501 centre hbration is indicated by the short arc straddhng the secondary (actu- 
ally a closed loop of width ~ 10"'^ in a^) while the cartesian motion along 
the epicycle path appears as a loop. The inset shows the evolution of the 

604 relative longitude A^, for 10^ yr including the period of libration around 0°. 

606 The example on the right panel shows the motion of MBA 121118 in a frame 

606 rotating with the mean motion of Vesta ('V') during a period of QS libration 
around Vesta lasting for 7 x 10^ yr. As in the previous case, the guiding cen- 
tre libration is indicated by the short arc straddling the secondary. Here, the 
higher amplitude of the fast harmonics in the evolution of act to smear 

610 out the epicycle to some extent. It is also evident in these plots that QS 

611 librators are physically located well outside the secondary's Hill sphere and 

612 should not be confused with keplerian satellites. 

613 Statistically, we find that 39 (24) out of the 129 (94) Ceres (Vesta) co- 

614 orbitals exhibited QS motion at some point during the 2 x 10^ yr period 
616 covered by the simulations reported in Section 3, a fraction of 25-30% in 

616 both cases. In the case of Vesta we find three episodes of QS libration of 

617 unusually long (> 4 x 10^ yr) duration. One is that of 121118 illustrated in 

618 Fig. [161 the others concern MBAs 22668 and 134633. In Fig. [T71 we compare 

619 the distribution of the relative proper semimajor axes of all objects that be- 

620 came temporary QS librators with the members of the persistent coorbital 

621 population. The two are generally separate with the former population fur- 

622 ther away and on either side of the latter. We believe this is because MBAs 

623 capable of becoming QSs are highly energetic. In other words, and referring 

624 to the top left and top right panels of Fig. [9l the value of the energy inte- 
626 gral 3C/8/i - represented by the horizontal line - is generally well separated 

626 from the potential S. Hence the amplitude of libration is only small near 

627 \r = 0° i.e. if the object is in QS mode at t = 0. This is the case for the 

628 one current and the two imminent QS orbiters of Ceres: 76146, 138245 and 

629 166529. The values of a,, for these objects are —0.0068, —0.0155 and —0.0083 

630 respectively. 

631 A particularly interesting subtype of QS librator is that with vanishing 

632 g uiding centre am plitude. These are the "retrograde satellites" (RS orbits) 



633 of lJacksonI (119131 ). In that case, the motion can be studied through an energy 

634 integral analogous to Eq. [T]but valid only in Hill's approximation to motion 
636 near the secondary. It depends explictly on e^, 7^, Ur and Qr (cf Eqs. 28 and 
636 29 ofN99). 
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Activating the QS mode requires ^ e while physical proximity of the 
object to the secondary in an RS orbit is controlled by e^, Ir and Ur- For 
the small values of e considered here, the relative eccentricity can still be 
small in absolute terms (e g lOe or 10~ ^). If Ir is also small (typically 



< Cr for bound orbits; see iNamounil . 119991 ). then such objects can remain, 
in principle, within a few times 10^ km of Ceres or Vesta. Note that small 
Cr implies a small libration amplitude for since < (see Section 3.3 
of N99). The long-term stability of these configurations depends, through 
the energy integral, on the evolution of the relative orbital elements. If these 
vary slowly, the asteroid remains trapped in QS motion for many libration 
cycles. 

To model the secular evolution of the asteroid's orbit we need a theoret- 
ical model of co-orbit al motion with in an N-body system. Such models do 
exist (eg the theory oflMessagd. Il966l valid near L4 and L5 and the more gen- 
eral model of iMoraisl . Il999l . l200ll ). However, the case in hand violates several 
assumptions for which those theories are strictly valid: low to moderate e 
and / of the asteroid (here e, / ^ e), a coorbital mode (QS) that was not ex- 
amined in those works and a clear separation of timescales between coorbital 
motion and the secular evolution of the orbit with w, Q <^ X (in fact, here we 
observe that X < w, Q). To check that the secular evolution of e, /, vu and Q 
of the asteroids does not depend on the presence of absence of the secondary 
mass we integrated the orbit of Vesta coorbital MBA 139168 with the eight 
major planets both with and without the massive asteroids (including Vesta) 
for 1 Myr. We found that the differences in the eccentricity and inclination 
vectors - eexp izu and /exp iQ respectively - are < 5% between the two cases. 
In addition, the asteroids' elements do not exhibit any of the features that 
arise from the three-body dynamics ( INamounil . Il999l ) . Hence, it is reasonable 
to assume that the secular forcing of and J,, is fully decoupled from the 
coorbital dynamics and can be modelled by N-body secular theory. 

In the absence of mean motion resonances, the secular evolution of the 
eccentricity and inclination vectors within a system of N bodies in near- 
circular, near-planar orbits around a central mass can be approximated by 
the so-called Laplace-Lagrange system of 2N first order coupled differential 
equations which are linear in e and / (Murray & Dermott 1999). The cor- 
responding secular solution for the eccentricity and inclination vectors of a 
particle introduced into that system has the form of a sum of N+l-periodic 
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573 complex functions: 



N 



686 



687 
688 



e = efexpi{gft + (3f) + ^Eiexpi{git + (3i) (7) 

i=l 
N 

I = J/expi (s/t + 7/) + ^/iexpi (sjt + 7i) (8) 



1=1 



574 where the first and second terms on the right-hand-side are referred to as 

676 the free and forced components respectively. The parameters of the forced 

676 component are derived directly from the Laplace-Lagrange solution. They, 

677 as well as the free eigenfrequencies gj and Sf, are dependent only on the 

678 semimajor axes of the planets and the particle. 

579 As with the theory of Morals, the validity of these expressions is limited to 

580 low-to-moderate e, / and gf ^ gi, Sf ^ Sj. Brouwer & van Woerkom (1950) 

581 showed that the 5:2 near-resonance between Jupiter and Saturn modifies 

582 the Laplace-Lagrange parameters of the real solar system. However, the 

583 independence of the forced component on the particle's e and I still holds; 

584 this allows us to write 

= e/expi(^/t + /3/) -e}expi(^}t + /3;^) (9) 
= J/expi(s/t + 7/) - J^expi + 7]:) (10) 



where the superscript refers to Ceres or Vesta as appropriate. Hence, 



686 and Ir are given by 



6-r — 6 4 



4 + ef + 2efe} cos [{gj - g}) t + {(3f - (3})] (11) 
Ir = lj + lf + 2IfI}cos[{sf-s})t+{^j-^})]. (12) 

The requirement for slow-evolving and Ir implies gf c:^ g'jr, Sf c:^ sj. Fur- 
ther, a necessary but not sufficient condition for and Ir to be small is 

689 Cf ^ ej, If ~ Ij. For these to be concurrently small, the following condition 

690 must also hold 

^ - A/3/ ^ Agf_ 

TT - A7/ ASf ^ ' 

591 where A/3j = /?/ - A7/ = 7/ - 7/, ^gf = gf - g} and As/ = s/ - s}. 

692 All these criteria, except the last one, can be tested for by making use of 

693 proper elements. The last criterion involves the proper phases which are not 
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594 included in the database but can, in principle, be recovered through harmonic 

595 analysis of the orbital element time series. 

596 As an example, we specify an upper limit in Ae/ and A sin // of 0.01 and 

597 a limit in Agj and Asf of 1 arscec per year. For Ceres co-orbitals, these 

598 correspond to a maximum epicycle excursion of ~ 0.05 AU and a period 

599 of > 1 Myr in the evolution of and J^. 15 and 7 of the 39 QS librators 

600 identified in this work satisfy either of these criteria respectively but none 

601 do both. On the other hand, 2 coorbitals of Vesta do satisfy both criteria (3 

602 and 8 respectively satisfy either one). Although neither of these two objects 

603 (78994 and 139168) satisfies the condition for concurrently small Cr and Ir 

604 we find the dynamical evolution of 139168 particularly interesting and we 

605 discuss it in some detail below. 

606 This asteroid has Ae/ = -0.0029379, A sin// = -0.0006901, Agf = 

607 0.020303 arcsec/yr and Asf = 0.075063 arcsec/yr. Inspection of our simu- 

608 lations of the nominal orbit of the asteroid shows that currently Cr ~ 0.013, 

609 sin/j, ~ 0.15. Both elements are slowly increasing in time from the present 

610 with a period significantly longer than the 2 x 10^ yr spanned by our inte- 

611 grations. Towards the past. The relative inclination continues to decrease 

612 with sinJr ~ 0.10 at t = —10^ yr while reaches a minimum of < 0.008 

613 at t = —8.5 X 10^ yr. For most of the integrated timespan, the object is a 

614 horseshoe with a very small opening angle, suggesting that it is sufficiently 

615 energetic to enter a QS mode (top left panel of Fig. ^ with 3C/8fi ~ 3.5. 

616 Indeed, we find two instances, indicated by the vertical arrows in the top 

617 panel of Fig. [181 when the asteroid is trapped into a QS mode for 10^ yr 

618 and an amplitude of ~ 1° (~ 10~^ rad). Since this is comparable to Cr 

619 (middle panel), these are not, strictly speaking, RS orbits. Nevertheless, it 

620 implies that the planar component of the motion occurs within a few times 

621 A^a of Vesta. In Fig. [12] we show the asteroid's motion for the phase of QS 

622 libration at t = —3.15 x 10^ yr in a heliocentric cartesian frame co-rotating 

623 with Vesta's mean motion. In the XY plane (bottom panel), the distance 

624 from Vesta varies between 10~^ and 8 x 10^^ AU. However, the excursion 

625 in Z (bottom panel) is significantly larger, ~ 0.3 AU. Note that the centre 

626 of the planar motion is offset from Vesta's position. This is partly because 

627 one is looking at the superposition of two harmonic modes (guiding centre 

628 and epicycle) of comparable amplitude. In addition, the potential minimum 

629 may not be exactly at A^ = since the maxima that bracket it are generally 

630 not equal unless Ur = kn. This becomes apparent if one expands the Hill 

631 potential to order higher than 2 or utilises the full potential of the averaged 
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632 motion i.e. Eq. [T] 

633 Since < Ir the QS-enabling potential minimum at the origin exists only 

634 for values of Ur sufficiently far from kir while the maxima on either side are 

635 highest for Ur = kir + jl^ becoming singularities if Ir = 0. These are the 

636 K > 2 orbits of N99. We can verify that this is the case here by overplotting 

637 Ur on the top panel of Fig. [TS] (dashed curve). We find that, during QS 

638 libration, the value of w,. is near ±90° (horizontal dotted lines) as expected. 

639 Eventually, and Ir will pass through their minima sufficiently closely 

640 in time to make long-term capture in a low amplitude QS mode possible. 

641 For the given values of Agf and Asf this should occur every ~ 6.4 x 10^ 

642 yr. A crude estimate on the expected number of such objects may be made 

643 under the assumptions that (a) this is the only object of its type, (b) the 

644 proper element catalog is complete for MBAs with a < 2.8 AU down to an 

645 absolute magnitude if ^ 16 and (c) the duration of QS capture is > 10^ 

646 yr, similar to large amplitude QS phases observed for other asteroids. The 

647 resulting average frequency of such objects at any one time is ^ 1.5 x 10~^ 



649 



648 Adopting an absolute magnitude distribution law of lO'^-''^ fiGladman et al. 



20091 ). we find that the frequency reaches unity for if ~ 26 or D = 12 — 37 

650 m objects. Unlike large amplitude Quasi-Satellites, the cartesian velocity of 

651 such objects with respect to Vesta is not high. For 139168 it ranges from 

652 2.2 to 3.2 km sec~^ for the integration spanning the last 1 Myr and will 

653 be lower if is smaller. Hence, and in view of possible in situ satellite 



654 searches by missions such as DAWN (jCellino et al.l . l2006l ) , newly discovered 

655 objects in apparent proximity to Vesta in the sky would need to be carefully 

656 followed up to determine whether they are "true" (ie keplerian) satellites or 

657 small-amplitude Quasi-Satellites. 



658 7. Conclusions and Discussion 

659 In this work we have demonstrated the existence of a population of Main 

660 Belt Asteroids (MBAs) in the coorbital resonance with the large asteroid 

661 (4) Vesta and the dwarf planet (1) Ceres. Libration within the resonance 

662 is transient in nature; our integrations show that these episodes can last for 

663 > 2 X 10^ yr. Partly due to the significant eccentricities and inclinations of 

664 these asteroids, we find that their dynamics are similar to those that govern 

665 the evolution of near-Earth asteroids in the 1:1 resonance with the Earth and 

666 Venus. However, due to the high density of objects as a function of semi- 

667 major axis, a steady state population of ~ 50 co-orb itals of Ceres and ~ 45 
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668 of Vesta is maintained. Apart from the natural dynamics of eccentric and 

669 inclined 1:1 librators, we identify two other mechanisms which contribute to 

670 the temporary nature of these objects. One is the inherent chaoticity of the 

671 orbits; with some exceptions, particles started in neighbouring orbits evolve 

672 apart over 10^ yr timescales. The other is close encounters with massive 

673 asteroids that do not participate in the Sun-Secondary-Particle three body 

674 problem. We show individual cases of asteroids leaving the coorbital reso- 
676 nance with Vesta following a deep encounter with Ceres. It is not clear if 

676 the latter effect is significant at the population level. It may be adversely 

677 affecting the occurence of long-lived tadpole librators of Vesta. Finally, we 

678 show that bound or Quasi-Satellite orbits around both Ceres and Vesta can 

679 exist and identify 3 current Quasi-Satellite librators around Ceres as well as 

680 an object which may experience episodes of long-lived bound motion within 
a few times 10~^ AU from Vesta. 

The demonstration of a resonant mechanism within the asteroid belt 
which acts independently of the major planets raises some interesting ques- 

684 tions to be addressed by future work. One of these is whether there exists a 

685 threshold below which a mass becomes a particle in real planetary systems. 

686 In partial response to this question - and as part of the integrations reported 

687 in Section 3 - we simulated the motion of several MBAs with proper semi- 

688 major axes within the coorbital region of 2 Pallas and 10 Hygiea. None of 

689 these was trapped in co-orbital libration which leads us to conclude that this 

690 threshold for the solar system's Main Asteroid Belt is ~ 10~^° solar masses. 

691 However, this may not be true of other planetary systems with fewer and/or 

692 less massive planets. We speculate that co-orbital trapping of planetesimals 

693 by terrestrial planets or large asteroids may give rise t o observationally verifi- 
able dynamical structures. This is not a new idea (eg Moldovan et al. . 2010l ) 



681 
682 
683 



694 



695 but our results show that even small bodies (10~^*^ < /i < 10~^) can maintain 

696 transient populations of Trojans in a steady state and that the dynamical 

697 excitation of planetesimals in a disk does not necessarily imply that the coor- 

698 bital resonance becomes ineffective. In this context, it may be relevant to the 

699 dynamical evolution of the larger planetesimals in protopl anetary disks and 



700 the p lanets or protoplanet cores embedded within them (jPapaloizou et al. 
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Motion of 139168 with respect to Vesta in cartesian ecliptic 
coordinates during an episode of QS capture at t ~ —3 x 10^ 
yr indicated in Fig. [181 The dashed hues mark the location 
X = Op, Vesta — 2.3615 au, y = z = 0. Note the large amplitude 
of the vertical motion ( right panel) compared to the planar 
projection (left panel) 
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Table 1: Dynamical parameters of large asteroids considered in this work. See text for 
details. 



Ceres Vesta Pallas Hygiea 

/i X 10^" 4.699 1.358 1.026 0.454 

e X 10^ 5.42 3.56 3.25 2.47 

Oproper (AU) 2.7670962 2.3615126 2.7709176 3.1417827 

eap™per(xlO^AU) 15.0 8.4 9.0 7.8 
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Table 2: Statistical results of the numerical simulations reported in Section 3. See text 
for details. 
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Table 3: Estimated mean and Full Width at Half Maximum (FWHM) of distributions 
presented in Fig. [51 
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Table 4: Results of the numerical simulations for the models described in Section 5. See 
text for details. 
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Figure 1: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 2: Cliristou and Wiegert 2010, Coorbitals of Ceres and Vesta 
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Figure 3: Cliristou and Wiegert 2010, Coorbitals of Ceres and Vesta 
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Figure 4: Christou and Wiegert 
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Figure 5: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 6: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 7: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 8: Christou and Wiegert 2010, Coorbitals of Ceres and Vesta 
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Figure 9: Christou and Wiegert 2010, Coorbitals of Ceres and Vesta 
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Figure 11: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 12: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 15: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 16: Christou and Wiegert 2010, Coorbitals of Ceres and Vesta 
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Figure 17: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 18: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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Figure 19: Christou and Wicgcrt 2010, Coorbitals of Ceres and Vesta 
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